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ABSTRACT 

This paper reports the search for quiescent low-mass X-ray binaries (qLMXBs) in the globular 
cluster (GC) NGC 6553 using an XMM-Newton observation designed specifically for that purpose. 
We spectrally identify one candidate qLMXB in the core of the cluster, based on the consistency of the 
spectrum with a neutron star H-atmosphere model at the distance of NGC 6553. Specifically, the best- 
fit radius found using the three XMM European Photon Imaging Camera spectra is i?NS = 6-3io'§ ^ m 
(for Mns = 1.4 M Q ) and the best-fit temperature is kT c g = 136^34 cV. Both physical parameters 
are in accordance with typical values of previously identified qLMXBs in GC and in the field, i.e., 
i?NS ~ 5-20 km and kT e g ~ 50-150 eV. A power-law (PL) component with a photon index T = 2.1+Qg 
is also required for the spectral fit and contributes ~ 33% of the total flux of the X-ray source. A 
detailed analysis supports the hypothesis that the PL component originates from nearby sources in 
the core, unresolved with XMM. The analysis of an archived Chandra observation provides marginal 
additional support to the stated hypothesis. Finally, a catalog of all the sources detected within the 
XMM field of view is presented here. 

Subject headings: stars: neutron — X-rays: binaries — globular clusters: individual (NGC 6553) 



1. INTRODUCTION 

The faint X-ray emission (L x ~ 10 32 -10 33 erg s _1 ) of 
transiently accreting low-mass X-ray binaries in quies- 
cent state (qLMXBs) was historically interpreted as a 
thermal blackbody with emission area smaller than that 
expected for a 10 km neutron star (NS). It was later 
claimed that the observed luminosity was not due to low 
accretion rate onto the N S surface, as initially suggested 
(jvan Paradijs et afll!987t ). but due to heat from the NS 
crust r adiating through th e upper layers of the NS atmo- 
sphere (jBrown et al.ll 19981) . In this interpretation, called 
Deep Crustal Heating, the energy is deposited in the 
NS crust by pressure-sensitive nuclear reactions (electron 
captures, neutron e missions and pycnonucle ar reactions, 
IGupta et all 120071: iHaensel fc Zdunikl [2001) . as matter 
accumulates on the surface during episodes of rapid ac- 
cretion. This chain of reactions, from the NS surface 
to depths with the density of undifferentiated equilib- 
rium, releases ~ 1.5 MeV per accreted nucleon in the 
crust. In a steady state, these reactions give rise to a 
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time-averaged luminosity, which is directly proportional 
to th e time-averaged mass accretion rate (jBrown et al.l 
1998): 



(L) = 9x10 



32 



(M) 



Q 



io- 



M Q yr" 1 1.5MeV/amu 



erg s 



(1) 

where Q is the average heat deposited in the NS crust 
per accreted nucleon. 

The resulting thermal spectrum is described by a re- 
alistic NS atmosphere composed exclusively of hydro- 
gen. Indeed, the gravitational settling of the accreted 
material from the low-mas s companion star ha ppens on 
time scales of —seconds (jBildsten et al.l I1992I ). result- 
ing in a pure H-atmosphere around the NS. The cur- 
rent models of NS H-atmosphere (Zavlin et al.l 119961 : 
iMcClintock etaTJ 12001 iHeinke et all 12006a!) show that 
the observed emission is consistent wit h that from the 
entire surface area of a -10 km NS (jRutledge et al.l 
1999). The H-atmosphere interpretation allows ob- 
servers to determine the physical radius i?NS (or the 
projected radius R x = i?Ns(l + z), with 1 + z = 

(l — 2GMns/-Rnsc 2 ) ) using spectral fitting in the 
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soft X-ray range, where the thermal spectrum of 
qLMXBs peaks. 

The radius measurement of NSs using the thermal 
spectrum of qLMXBs is one way to place constraints 
on the equation of state (EoS) of dense nuclear matter. 
This relation between the pressure and the density in the 
interior of NSs is highly uncertain and multiple compet- 
ing theories are proposed to explain the behavior of cold 
matter at densities above 2.35x 10 14 g cm' 3 . Quiescent 
LMXBs found in globular clusters (GCs) are routinely 
used to produce precise measurements of NS radii. The 
distance to GCs are generally known with precisions of 
the order of ~ 5%-10% and the overab undances of NS 
binary systems in GCs ( Hut et al.lll992l ) make these ob- 
jects ideal targets to pursue searches of qLMXBs capable 
of providing useful constraints on the dense matter EoS. 

Twenty-six GC qLMXBs, including the more recent 
candidates, have been discov ered so far (see lHeinke et al.l 
I2003al iGuillot et all l2009ai for complementary lists). 
However, for some of them, the poorly constrained _Rns 
and kT e g measurements, as well as the high galac- 
tic absorption in their direction (for example, -/Vh ~ 
1.2 xlO 22 atoms cm 2 in the direction of Terzan 5), af- 
fects the certainty of their identification. While qLMXBs 
in the field of the Galaxy or in GC have been observed in 
quiescence following outbursts events in historical X-ray 
transients, none of the spectrally identified qLMXBs in 
GC (with spectra consistent with NS atmosphere mod- 
els or via their X-ray spectral colors) have been ob- 
served in outbursts. There has been one exception re- 
cently discovered. An X-ray transient i n the cluster 
Terzan 5 has been obser v ed in outburst ([Bordas et al.1 
[20101 iPoolev et all [2010l : iDegenaar fc Wiinandsl l20TTf T 
after being tentativ ely classified as a qLMXB (CX25, 
IHeinke et al.ll2006blh To increase the list of known GC 
qLMXBs, a program of short observations using the Eu- 
ropea n Photon Imaging Cameras (EPIC; iStruder et al.l 
120011 ) onboard XMM-Newton have been undertaken to 
survey GCs and search for qLMXBs. This paper reports 
the discovery via spectral identification of a candidate 
qLMXB in the core of NGC 6553. A short archived 
Chandra observation targeted at this GC was also ana- 
lyzed to provide more accurate positions and tentatively 
confirm the source classification. 

The targeted cluster, NGC 6553, is a low-galactic 
latitude GC located at the position R.A.=18 h 09 m 17.6 s 
and deck = -25 54'31'/3 (J2000), approximately 2.2 kpc 
from the galactic center, corresponding to a heliocen- 
tric distance of d = 6.0 kpc. It is a GC of moder- 
ate core compactness and moderate size: core radius 
r c = 0.'55, half mass radius thm = 1.'55 an d tidal radius 
r t = 8.16. Its metallicity [M/H] = -0.09 (jValenti et al.1 
2007) makes it one of the most metallic clusters in the 
Galaxy. The foreground reddening E (B — V) = 0.63 
corresponds to a moderately high hy drogen column den- 
sity A H = 0.35 xlO 22 atoms cm -2 (jPredehl fc Schmittl 
1995). The hydrogen column density is written ATh,22 = 
0.35 afterward, and this value will be used for the X-ray 
spectral fits. This value is similar to th e value found from 
NRAO data (jDickev fc Lockmanlll990D 1 , which accounts 
for the atomic hydrogen Ahi through the entire galaxy. 
The GC properties come from the catalog of GCs (jHarris! 

1 From http://cxc.harvard.edu/toolkit/colden.jsp 



119961 update Dec. 2010), except when other references 
are provided. 

The organization of this paper is as follows. Section [2] 
describes the data reduction and analysis of the XMM- 
Newton and Chandra observations. In Section [3l we 
present the results and we discuss them and conclude 
in Section @] 

2. DATA REDUCTION AND ANALYSIS 

NGC 6553 was observed with XMM-Newton on 2006 
October 6 at 01:41:44 UT, using the three EPIC cameras, 
for an exposure time of 20.4 ks, with the medium filter. 
The Chandra X-ray Observatory observed the GC for 
5.25 ks, on 2008 October 30 at 01:38 :12 UT, with the 
ACIS-S detector (|Garmire et al.ll200l in VERY-FAINT 
mode. 

2.1. Data reduction 

Reduction of XMM-Newton data is performed with the 
XMM-Newton Science Analysis Software (SAS) v8.0.0, 
using the standard procedures 2 . More specifically, 
epchain and emchain are used for the preliminary re- 
duction of the raw pn and MOS data files, respectively. 
Single and double patterns are used for pn data and sin- 
gle, double and quadruple patterns for the MOS cam- 
eras, both in the 0.3-10 keV range. Prior to the source 
detection, the data are checked for background flares by 
looking for > 3<r deviation from the mean count rate of 
the entire detector. None are found and the whole inte- 
gration time is used. Also, an exposure map is created to 
correct for the effe ct of vignetting. The s cript wavdetect 
from CI AO v4.2 (jFruscione et al.ll2006fi is then run on 
the pn image in the 0.3-10 keV range, with the follow- 
ing parameters: minimum relative exposure of 0.1, the 
wavelet scales "1.0 2.0 4-0 8.0" and a significance thresh- 
old 3xl0 -6 . The latter is equal to the inverse of the 
number of pixels in the image, and allows for about one 
spurious detection in the entire image. Sources with de- 
tection significance a > 4 are retained for analysis. 

The reduction and analysis of the Cha ndra data are ac- 
comp lished with the CI AO v4.2 package (jFruscione et al.1 
2006). The level- 1 event file is first reprocessed using 
the public script chandra_repro which performs the steps 
recommended by the data preparation analysis thread 3 
(corrections for charge transfer inefficiency, destreaking, 
bad pixel removal, etc, if needed) making use of the latest 
effective area maps, q uantum efficiency ma ps, and gain 
maps of CALDB v4.3 (jGraessle et al.ll2007t) . The newly 
created level-2 event file is then checked for background 
flares, but there are none detected. 

2.2. Count extraction 

The command evselect from the SAS is used for 
source and background count extraction of the XMM 
data. A radius of 30" around each source is chosen, 
accounting for 88% of the total energy from an on-axis 
source at 1.5 keV with the pn detector 4 . The background 
counts, used for background subtraction, are extracted 
from a larger region (100") around the source of interest, 

2 User Guide to the XMM-Newton Science Analysis System, 
Issue 5.0, 2008 (ESA: XMM-Newton SOC). 

3 http : //exc .harvard. edu/ciao/threads/data. html 

4 From XMM-Newton Users Handbook, Figure 3.7, July 2010 
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excluding the source itself and other detected sources in 
close proximity with a circular region of 40" around each 
of them, ensuring that 92% of the overlapping source 
counts (at 1.5 keV on the pn detector) are excluded. In 
some cases, the close proximity of the detected sources 
requires to adapt the exclusion radii of nearby sources. 
Those few cases are discussed in Section [31 The tasks 
rmfgen and arfgen generate the response matrix files 
(RMF) and the ancillary response files (ARF) for each 
observations. The XMM-Newton extracted counts are 
binned using grppha into energy bins between 0.3 and 
10 keV with a minimum of 20 counts per bin. This cri- 
terion ensures having approximate Gaussian statistics in 
each spectral bin. 

For Chandra data, the script psextract is used to ex- 
tract the counts of the X-ray sources, in the energy range 
0.5-8 keV. The extraction region depends on the off-axis 
angle since the point spread function (PSF) degrades at 
large off-axis angles. Since the analyzed observation of 
NGC 6553 was performed with the ACIS-S instrument 
and a focal plane temperature of -120 °C, the RMFs have 
to be recalculated, according to the recommendations of 
the CIAO Science Thread "Creating ACIS RMFs with 
mkacisrmf " . It is also crucial to recalculate the ARFs 
using the new RMFs in order to match the energy grids 
between the two files. 

In both the XMM-Newton and Chandra observations, 
the count rates are low enough for the effects of pile-up to 
be ignored. Specifically, for the Chandra data, the count 
rate was 0.031 photons per detector frame (3.24 sec), and 
therefore the effect of pile-up can be safely neglected. 

2.3. Spectral Analysis 

The identification of candidate qLMXBs in this analy- 
sis is based on the spectral identification with the mod- 
els that empirically describe the spectra of previously 
known qLMXBs. Spectra from the pn camera are first 
extracted and analyzed, since its superior efficiency facil- 
itates spectral analysis with an improved signal-to-noise 
ratio (S/N) over the two MOS c ameras. H-atm osphere 
models are available in XSPEC (1 Arnaudl 1 1 9961 ) . includ- 
ing nsa and nsagrav ([Zavlin et al.l 119961) an d nsatmos 
(|McClintock et all 12004 IHeinke et al.ll2006al) . The ma- 
jor difference between the models resides in the sur- 
face gravity values g used to calculate the model: 
nsa has been implemented with a fixed value g — 
2.43xl0 14 cms -2 , nsagrav has been computed for a 
range g = (0.1-10) x 10 cms" 2 and nsatmos for the 
range g = (0.63-6.3) xlO 14 cms -2 . Previous works 
have demonstrated that NS atmosphere models calcu- 
lated with multiple values of g ar e better adapted to 
fit the thermal spectr a of qLMXBs (Heinke et aL 2006a; 
iWebb fc Barretll2007D . We choose here to use nsatmos. 

The identification of candidate qLMXBs is based on 
the statistical consistency (null hypothesis probability 
n.h.p> 10~ 2 ) with the nsatmos model at the distance of 
the host GC, with the galactic absorption taken into ac- 
count using the multiplicative model wabs. Additional 
conditions impose that the best-fit effective tempera- 
ture and projected radius are within the range of previ- 
ously observed GC qL MXBs, i.e., kT eS ~ 5 0-180 eV and 
floo ~ 5-20 km (see IHeinke et all I2003bt iGuillot et all 
2009a, for complementary lists of GC qLMXBs). For 



identified candidate qLMXBs, the pn, MOS1 and MOS2 
extracted spectra are then fit simultaneously to improve 
statistics and diminish the uncertainties on the best-fit 
parameters. 

In some cases, an excess of counts at high energy 
is modeled with an additional power-law (PL) compo- 
nent. This is empirically justified by the fact that 
qLMXBs sometimes display a hard PL which domi- 
nates the spectrum above 2 keV. Most qLMXBs in 
the field of the galaxy have spectra best fitted with H- 
atmosphere models co mbined with a hard-photon PL 
(for ex ample, Ce n X-4. lAsa et aL 119961: [k utlcdge et all 
2001U Aql X-l lAsai et all 119981: iCampana et all 119981: 



Rutledge et alll2001a|) . ithas been claimed that some 



GC qLMXBs also displayed significant PL contribu- 
tions (between ~ 10% and ~ 50% of the total flux), 
for examp le, in NGC 6440 (ICackett et "all l2005h. and 
in Ter 5 (Heinke et all l2006bl ) . Proposed interpreta- 
tions of the observed PL tail in clude residual accre- 
tion onto the NS m agnetosphere (jGrindlav et all 120011 : 
ICackett et alll2005l ), shock emission via the emergence 
of a magnetic field (jCampana fc Stella 2000), or an in- 
trabinary shock between the winds fro m the NS and its 
companion star (jCampana et al.ll2004l ). However, analy- 
ses of the quiescent emission of LMXBs have shown evi- 
dence of the presence of a variabl e low-level accretion o n 
the NS; for the LMX B Aql X-l (iRutledge et all 120021), 
for XTE J1701-462 dFridriksson et allfcoiOft and for 



the LMXB Cen X-4 (jCackett et al.ll201oh . 

Such an additional spectral component for qLMXBs 
in the core of GCs might not necessarily originate from 
one of the interpretations listed above. For example, 
a PL component was required in addition to the H- 
atmosphere model for the spectral fit of the candidate 
qLMXB XM MU J171433-292747 in NGC 6304 observed 
with XMM (jGuillot et all l2009al ). However, a short 
5ks archived Chandra observation demonstrated that 
the high energy excess of counts noticed for this source 
originated from another close by X-ray source (possibly 
a cataclysmic variab le) unresolved in the XMM data 
(|Guillot et al.ll2009bt ). Then, a longer 100 ks Chandra 
exposure confirmed the NS radius measurement and the 
initially reported results (S. Guillot et al. 2011, in prepa- 
ration) . 

The main focus of this paper is the discovery of candi- 
date qLMXBs in NGC 6553, but as part of the systematic 
analysis of XMM data, a simple spectral analysis is per- 
formed for all other X-ray sources that are not consistent 
with NS H-atmosphere models at the distance of the host 
GC. Specifically, a simple absorbed PL model is used and 
their best-fit photon indices, as well as their unabsorbed 
flux (0.5-10 keV) are reported. In some cases, the X-ray 
spectra cannot be fit adequately with either an NS H- 
atmosphere model or a PL model. In those cases, the 
use of other models is attempted for completeness and 
the details are provided in Section [3] Finally, sources 
with an S/N < 3 are fit with a fixed (r = 1.5) photon 
index. 

In all cases, the galactic absorption is taken into ac- 
count with the model wabs, keeping the value of 2Vh fixed 
at the value in the direction of the GC, i.e., Ah, 22 = 0.35. 
Also, all quoted uncertainties on the spectral parameters 
are 90% confidence level in the text and the tables. Fi- 
nally, the quoted fluxes are corrected for the finite aper- 
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Fig. 1 — This EPIC/pn image of NGC 6553 shows the 43 sources 
detected on the pn camera. The small and large dashed circles 
respectively represent the half-mass (rjjM = 1/55) and tidal radii 
(rt = 8.' 16). A log gray-scale has been used for the image due to 
the brightness of source #1. In this figure, north is up and R.A. 
increases to the left. The candidate qLMXB detected in this GC 
is the source #3, corresponding to XMMU J180916-255425. 

ture (extraction radius), except noted otherwise. 

2.4. Astrometric correction 

The uncertainties reported in Table Q] on the XMM 
positions are composed of the wavdetect statistical 
error on the pn positions (up to 1") and the systematic 
uncertainty of XMM (2"at lc 5 ). An absolute astromet- 
ric correction - into the frame of the Two Micron All Sky 
Survey (2MASS) - is applied through the association of 
X-ray sources with 2MASS counterparts. Three X-ray 
sources were firmly associated with 2MASS objects. 
The sources #5, #2 and #9 have probabilities 0.001%, 
0.21% and 0.04%, respectively, that a star as bright or 
brighter as the stars 2MASS J18093564-25555377, 
2MASS J18100723-25511760 and 

2MASS J18083961-2601203, respectively, lie as close or 
closer to the three X-ray sources. These probabilities 
were calculated using all 2MASS stars detected in 2' 
wide annulli, centered at the GC optical center, and 
with respective average radii equal to the distances 
of the respective 2MASS counterparts from the GC 
center. This means that sources #5, #2 and #9 have 
probabilities 99.999%, 99.79% and 99.96% to be asso- 
ciated with the stars 2MASS J18093564-25555377, 
2MASS J18100723-25511760 and 

2MASS J18083961-2601203, respectively. From 
these three associations, the best-fit required transfor- 
mation is a shift of -O'.'l in R.A. and 0'.'71 in decl. The 
residual uncertainties are 0"34 in R.A. and 1'.'08 in decl. 
This is added in quadrature to the 2MASS astrometric 

5 f rom XMM Scien ce Operations Centre XMM-SOC-CAL-TN- 
0018. [GlTalmzzil 112008) 



uncertainty (~ 0'.'15 6 ) and to the statistical uncertainty 
of the X-ray source detection. 

3. RESULTS 

In the reduced image of NGC 6553, the wavdetect 
algorithm detected 43 X-ray sources on the pn camera 
within the field of the observation. Table [T] lists all the 
sources detected with their positions, their statistical un- 
certainty and their detection significance, while Table [2] 
provides the results of the spectral fits and the fluxes. 
The spectral analysis of candidate qLMXBs and of the 
sources for which a PL model does not describe the data 
is detailed below. 

As can be seen in Figure [IJ the proximity (17'.' 5) of 
source #15 to source #5 requires to adapt the extrac- 
tion radius prior to the spectral analysis. More specifi- 
cally, instead of excluding the default 40" around nearby 
sources, the exclusion radii of source #15 is chosen to 
be 10" (59% of encircled count fraction, ECF) for the 
count extraction of #5, and the exclusion radii of source 
#5 is chosen to be 12"5 (68% of ECF) for the count ex- 
traction of source #15. The radii mentioned were deter- 
mined by visual inspection of the X-ray image. Similarly, 
the wavdetect algorithm detected two sources, #11 and 
#14, located 16" apart near the edge of the pn camera. 
In this case, the exclusion radii of #11 and #14 are set 
to 15" and 10" 5, when considering the extraction region 
of sources #14 and #11, respectively. 

3.1. Source #1 - XMMU J 180846- 260043 

This X-ray source has the largest pn count rate of the 
observation, 0.197 ± 0.004 cts s~ x . Nevertheless, it is not 
affected pile-up (accumulation of photon events within a 
single CCD time frame, 73 ms) since its count rate per 
detector frame is then ~ 0.014. 

The bright X-ray source corr esponds to the RO SAT 
source 1RXS J180847.7-260037 (jVoees et al.l[l99l and 
to the IR star 2MASS J18084678-2600417 (with prob- 
ability of association 99.999%). It is also named, 
HD 315209, as listed in the Tycho Reference Catalog 
with the proper motion +14 .2±2.4 mas yr -1 in R.A . and 
-g.OilJmasyr- 1 in decl. (IHog et all 1 1991 l2000h . The 
star has a spectral type K5V and is located at a dis- 
tance of 33 pc, de duced from spectroscopic observations 
(jRiaz et al.ll200ll . 

The XMM source is adjacent to a bad column on the 
pn camera. Therefore, the resulting astrometry is likely 
to be offset in the direction opposite to the bad column, 
i.e., toward increasing right ascension. Therefore, we 
provide the position obtained with the MOS2 camera 
- R.A.=18 h 08 m 46.7 s and decl. = -26°00'42"7. 

The high-S/N spectrum of this star is not well 
described by an absorbed PL: x^/dof (prob.) = 
1.80/120 (1.7xl0~ 7 ). Removing the absorption, i.e., 
setting A r H=0, also leads to a non- acceptable fit with 
n.h.p. ~ 10~ 60 . Similarly, NS H-atmosphere, ther- 
mal bremsstrahlung and blackbody models do not fit 
properly the data, leading to n.h.p. < 10~ 10 , whether 
or not galactic absorption is included with the mod- 
els fitted. A Raymond-Smith (RS) plasma model pro- 
vides improved fit statistics with n.h.p. m 10 -3 but 

e Section 2 of http: //www. ipac . caltech.edu/2mass/releases/allsky/doc/ 
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TABLE 1 

X-RAY SOURCES DETECTED IN THE XMM OBSERVATION OF NGC 6553. 



Object Name 


RA 


decl. 


5R.A.\<Wi. a 


S/N b 


ID 








(J2000) 


(J2000) 


(") 








XMMU 


J180846- 


-260043 


272.19551 


-26 


.01220 


±0.0\0. 





166.1 


01 


XMMU 


J181007- 


-255118 


272.53028 


-25 


.85518 


±0.3\0. 


3 


28.0 


02 


XMMU 


J180916- 


-255425 


272.31889 


-25 


90704 


±0.4\0. 


3 


26.2 


03 


XMMU 


J180832- 


-255852 


272.13472 


-25 


.98135 


±0.4\0. 


1 


24.6 


04 


XMMU 


J180935- 


-255554 


272.39844 


-25 


93172 


±0.4\0. 


5 


22.6 


05 


XMMU 


J180845- 


-260035 


272.18903 


-26 


00978 


±0.4\0. 


•1 


21.4 


06 


XMMU 


J180935- 


-255157 


272.39592 


-25 


86605 


±0.3\0. 


3 


21.3 


07 


XMMU 


J180952- 


-254637 


272.46964 


-25 


.77697 


±0.7\0. 


6 


16.9 


08 


XMMU 


J180839- 


-260118 


272.16497 


-26 


02191 


±0.6\0. 


5 


16.7 


09 


XMMU 


J180936- 


-255952 


272.40357 


-25 


99803 


±0.6\0. 


5 


15.5 


10 


XMMU 


J180851- 


-254113 


272.21633 


-25 


.68706 


±0.6\0. 


5 


14.1 


11 


XMMU 


J180910- 


-254125 


272.29530 


-25 


69036 


±0.5\0. 


5 


12.1 


12 


XMMU 


J180931- 


-254105 


272.38234 


-25 


.68474 


±0.6\0. 


1 


10.9 


13 


XMMU 


J180853- 


-254106 


272.22092 


-25 


68513 


±0.6\0. 


5 


9.3 


14 


XMMU 


J180934- 


-255549 


272.39324 


-25 


93030 


±0.5\0. 


4 


9.2 


15 


XMMU 


J180915- 


-255631 


272.31424 


-25 


.94212 


±0.7\0. 


5 


9.0 


16 


XMMU 


J180921- 


-260325 


272.34126 


-26 


.05712 


±1.0\0. 


7 


8.0 


17 


XMMU 


J181002- 


-260135 


272.51131 


-26 


02663 


±0.6\0. 


7 


7.9 


18 


XMMU 


J180922- 


-255310 


272.34343 


-25 


. oouoo 


±0.9\0. 


6 


7.5 


19 


XMMU 


J180854- 


260353 


272.22729 


-26 


.06497 


±0.6\0. 


5 


7.4 


20 


XMMU 


J180937- 


-255408 


272.40664 


-25 


.90245 


±0.5\0. 


5 


7.3 


21 


XMMU 


J180912- 


-255138 


272.30407 


-25 




±0.8\0. 


8 


6.6 


22 


XMMU 


J180931- 


-255149 


272.38286 


-25 


.86367 


±0.8\0. 


6 


6.4 


23 


XMMU 


J180915- 


-260127 


272.31512 


-26 


.02427 


±0.9\0. 


9 


6.3 


24 


XMMU 


J180938- 


-260233 


272.41026 


-26 


04259 


±1.0\0. 


6 


6.2 


25 


XMMU 


J180856- 


-254806 


272.23482 


-25 


.80185 


±0.9\0. 


7 


6.2 


26 


XMMU 


J180921- 


-255731 


272.34030 


-25 


.95878 


±0.8\0. 


7 


6.2 


27 


XMMU 


J180920- 


-255614 


272.33388 


-25 


93725 


±0.7\0. 


6 


5.9 


28 


XMMU 


J180917- 


-254614 


272.32096 


-25 


.77083 


±0.9\0. 


7 


5.8 


29 


XMMU 


J180941- 


-255048 


272.42476 


-25 


.84671 


±0.9\0. 


8 


5.7 


30 


XMMU 


J180913- 


-254834 


272.30522 


-25 


.80972 


±0.9\0. 


8 


5.6 


31 


XMMU 


J180943- 


-255808 




-25 


.96905 


±0.7\0. 


7 


5.6 


32 


XMMU 


J180950- 


-255531 


272.45886 


-25 


.92548 


±0.7\0. 


6 


5.5 


33 


XMMU 


J180930- 


-260205 


272.37870 


-26 


.03489 


±0.8\0. 


8 


5.5 


34 


XMMU 


J180939- 


-254723 


272.41319 


-25 


.78992 


±0.0\0. 





5.4 


35 


XMMU 


J180950- 


260307 


272.46218 


-26 


.05219 


±0.9\0. 


7 


5.4 


36 


XMMU 


J180839- 


-254707 


272.16615 


-25 


.78536 


±0.7\0. 


6 


5.3 


37 


XMMU 


J180934- 


-254537 


272.39339 


-25 


.76053 


±0.6\0. 


6 


5.2 


38 


XMMU 


J181003- 


-255636 


272.51627 


-25 


.94347 


±0.7\0. 


6 


5.0 


39 


XMMU 


J180844- 


-254812 


272.18444 


-25 


.80336 


±0.8\0. 


7 


4.6 


10 


XMMU 


J180924- 


-255540 


272.35017 


-25 


.92790 


±0.7\0. 


6 


4.5 


41 


XMMU 


J181007- 


-260503 


272.53191 


-26 


.08422 


±0.7\1. 





1.1 


42 


XMMU 


J181023- 


-254811 


272.59968 


-25 


.80325 


±0.7\0. 


7 


4.2 


43 



a Statistical uncertainty on the position 

k Detection significance provided by wavdetect 

this is still not a statistically acceptable fit. A two- 
temperature RS model, often used t o represent stellar 
coronal spectra (|Dempsev et al.|[l993f) . satisfactorily de- 
scribes the spectrum, xt/doi (prob.) = 1.01/11.5 (0.44), 
with best-fit temperatures fcT c fF .1 = 0.35^0 x2 kcV and 
kT cSa = 2.80i^fkeV. With this model, the best-fit 
hydrogen column density is iVn,22 = 0.09 ± 0.04, consis- 
tent with a small distance to this star which lies in the 
foreground of NGC 6553. 

3.2. Source #5 - XMMU J180935-255554 

This source, located within the GC tidal radius, is 
likely associ ated with the star HP 165 973 of spectral 
type G8III (|Houk fc Smith-Moord H9881 The 2MASS 
name of this star is 2MASS J18093564-2555537 (99.999% 
chance of association with the X-ray source #5). It 
has a proper motion of +12. 9zfcl.9masyr~ 1 in R.A . and 
-S.Sil.emasyr^ 1 in decl. (Tffog et al.l " ll998l l2000h . No 
parallax measurement is available for this star. 



The X-ray spectrum cannot be fit with the NS H- 
atmosphere, PL, blackbody and thermal bremsstrahlung 
models. All of these models lead to n.h.p. < 0.006. How- 
ever, an RS plasma model provides a statistically accept- 
able fit with xl/doi (prob.) = 0.75/13 (0.71). The best- 
fit temperature is kT c g — O.S?^'^ keV and the best-fit 
hydrogen column density is -/Vh,22 < 0.26 (90% confi- 
dence upper limit). 

3.3. Source #9 - XMMU J180839-260118 

The source #9 has a soft thermal spectrum that is 
found to be acceptably fit by the nsatmos model: fcT er j = 
80ii2 e V and R NS = 15.4±^km for M NS = 1.4 M s , 
with x?/dof (prob.) = 1.50/32 (0.03). The spectrum is 
also acceptably fit with a RS plasma model: kT e g = 
0.87±°;^keV, iV H ,22 = 0.19^8 and Z = 0.10±°;^ 
with x^/dof (prob.) = 1.20/30 (0.21). A Kolmogorov- 
Smirnov test (probability p — 0.22) shows no evidence of 
variability. 
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TABLE 2 

Spectral results of X-ray sources in NGC 6553. 



ID 


pn Rate 


r 


kT cS 


Fx xS/d.o.f. (prob.) 


Model 




(cts ks" 1 ) 




(keV) 




Used 



01 


197.3±3.8 


_ 


0.35+H1 ; 2.80+^ 5 7 


17.3 


1.01/115 


(0.44) 


RS+RS 


02 


17.7±1.3 


3.0+2-3 




2.54 


1.32/15 


(0.18) 


PL 


03 


27.1±1.5 a 


9 1+0.5 

8 


o.i36tg;j»i 


1.29 


0.90/54 


(0.68) 


NSATMOS+PL 


04 


18.1±1.4 


I 4+0-2 
— 2 




2.67 


0.79/18 


(0.71) 


PL 


05 


15. Oil. 3 






0.42 


0.75/13 


(0.71) 


KS 


06 


8.6±0.9 


3 3+ - 5 

— 0.4 




5.89 


1.50/ 7 


(0.16) 


PL 


07 


17.6±1.4 


u ' a — 0.2 




2.46 


0.75/18 


(0.76) 


PL 


08 


10.8±1.1 


u -°-0.3 




2.72 


0.94/11 


0.50) 


PL 


09 


11.6±1.2 




n 873+°' 145 
U.O(O_ 092 


0.65 


1 20/30 


0.21) 


RS 


10 


9.6±1.1 


2 4+° ! ! 




0.74 


0.75/12 


0.70) 


PL 


\\ 


6.0i0.9 


n 7 +0.3 
u - '-0.4 




2.58 


0.85/ 6 


0.53) 


PL 


12 


4.0±1.0 


n o + 0.7 
-0.3+07 




2.44 


0.83/ 8 


0.58) 


PL 


13 


3.9±0.9 


i 9+O.6 
L -*-0.6 


- 


1.07 


0.92/ 6 


(0.48) 


PL 


11 


3.7±0.7 


n 9 +0.6 


- 


2.2 


0.66/ 3 


(0.57) 


PL 


15 


6.7±1.0 


9 q+0.8 


- 


0.44 


2.02/ 8 


(0.04) 


PL 


16 


4.5±1.0 


4-5±i;l 


- 


0.31 


1.27/ 9 


(0.25) 


PL 


17 


4.3±0.9 


3.9+1? 


— 


0.34 


2.15/ 7 


(0.03) 


PL 


18 


5.1±1.0 


1.5+2 : 5 

— 0.5 




0.96 


0.59/ 8 


(0.79) 


PL 


19 


3.7±1.0 


1.3+2-1 

— 0. * 




0.30 


0.86/ 9 


(0.56) 


PL 


20 


4.0±0.9 


5.0+11 

^ ^ —1.2 




0.51 


1.30/ 7 


(0.25) 


PL 


21 


4 4-1-0.9 


3 5+ 10 




0.31 


1.14/ 7 


0.33) 


PL 


22 


2.3±0.9 


(1.5) 




0.13 


0.54/ 8 


0.82) 


PL 


23 


4.6±1.0 


2.1+nI 
—0.6 




1.60 


0.65/ 8 


(0.74) 


PL 


9/1 


o.ozcu.y 


i n+°- 6 

i - U -0.6 




ft A7 


1 .01/ ( 


/ft OA\ 


i 


25 


4.4±0.9 


- .-1-17 


— 


0.44 


1.84/ 7 


(0.08) 


PL 


26 


2.6±0.9 


1 4+ 12 


- 


0.28 


0.52/ 6 


(0.79) 


PL 


27 


1.5±0.9 


(1.5) 




0.08 


0.95/ 8 


(0.37) 


PL 


28 


0.4±0.9 


(1.5) 




0.05 


2.26/ 7 


(0.03) 


PL 


29 


3.3±0.8 


n S +0.9 




0.75 


0.34/ 5 


(0.89) 


PL 


30 


3.0±0.8 


1 s+ 13 




0.29 


0.60/ 5 


(0.70) 


PL 


31 


3.0±0.9 


1 4+0-9 
l - q: -0.8 




0.30 


0.36/ 7 


(0.93) 


PL 


32 


1.5±0.9 


(1.5) 




0.20 


0.39/ 7 


(0.91) 


PL 


33 


2.2±0.9 


(1.5) 




0.21 


0.47/ 7 


(0.86) 


PL 


34 


2.7±0.8 


1 6+ 12 
J - u — l.i 




0.31 


0.56/ 5 


(0.73) 


PL 


35 


4.1±0.8 


n O + 0.7 

o.8_ Q 7 




0.51 


0.63/ 8 


(0.85) 


PL 


36 


3.8±0.9 


2.8i°; 7 




0.45 


1.66/ 7 


(0.11) 


PL 


37 


2.0±0.8 


(1.5) 




0.09 


2.25/ 6 


(0.04) 


PL 


38 


3.3±0.9 


3 o+2-l 




0.33 


1.07/ 6 


(0.38) 


PL 


39 


3.0±0.9 


4 3+3-4 




0.33 


0.63/ 6 


(0.70) 


PL 


10 


1.7±0.8 


(1.5) 




0.20 


0.66/ 5 


(0.65) 


PL 


41 


0.7±0.9 


(1.5) 




0.01 


2.01/ 7 


(0.05) 


PL 


12 


0.4±0.7 


(1.5) 




0.24 


1.37/ 4 


(0.24) 


PL 


13 


1.8±0.8 


(1.5) 




0.49 


0.88/ 5 


(0.50) 


PL 



Note. — The columns of this table are, from left to right: the source ID according to Tabled the 
background subtracted count rates in the 0.3—10 keV range, the best-fit photon index when a power law was 
used (values in parentheses were held fixed), the best-fit effective temperature when a thermal model was 
used, and the unabsorbed X-ray flux in unit of 10 — 13 erg cm -2 s — 1 in the band 0.5— lOkeV. Note that the 
flux is not corrected for the finite aperture in this table. This energy range is different from that used for 
the spectrum extraction (0.3-10 keV), for comparison purposes with previously published work. The last 
column is the resulting X^-st&t istic and the model used (PL: power law, RS: Raymond-Smith plasma). All 
fits have been performed with a fixed value of the galactic absorption, Nu,22— 0.35. 

a The MOSl/2 count rates are respectively S.5±0.8 and 6.6±0.7 counts per ksec in the 0.3-10 keV range. 
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Fig. 2. — XMM-pn spectrum of the candidate qLMXB 
XMMU J180916-255425 in NGC 6553. The two components (H- 
atmosphere nsatmos and PL powerlaw) arc shown in dashed lines. 
They are both affected by interstellar absorption modeled with 
wabs and the hydrogen column density Afjj,22=0.35. The PL dom- 
inates the spectrum at E}t 2keV. 

The infrared (IR) source, 2MASS J18083961-2601203, 
is located at a distance of l'/51 from source #9 (after 
correction) and is identified as the probable counterpart 
of the candidate qLMXB, with a probability of associa- 
tion of 99.86%. It has magnitudes in the J, H and K 
IR bands of mj = 10.146(15), m H = 9.796(23), and 
niK = 9.723 (31), where the numbers in parentheses in- 
dicate the uncertainties on the magnitudes. The X-ray to 
J-band flux ratio of the association is Fx/Fj ps 2xl0 -3 
and is consistent with that of a coronally active star. The 
possibility that source #9 is a qLMXB with a giant star 
in NGC 6553 cannot be discarded on this basis. 

However, in the col or-magnitude diagram of NGC 6553 
(jOrtolani et al.|[l995l ). the y- band magnitude of this star 
in the NOMAD Catalogue (|Zacharias et all I2005T ) lies 
~ 5 mag above the tip of the red giant branch. This 
observation strongly disfavors the giant star hypothesis, 
since the star is ~ 100 x brighter than the brightest gi- 
ants in the cluster. In addition, the location of the sys- 
tem outside the tidal radius of the cluster also disfavors 
the qLMXB classification. Therefore, we conclude that 
source #9 is not a candidate qLMXB at the distance of 
the GC. 

3.4. Candidate qLMXB, source #3 - 
XMMU J 180916- 255425 

3.4.1. XMM-Newton observation 

The X-ray spectrum of this source is acceptably fit 
(x^/dof (prob.) = 1.44/33 (0.05)) with nsatmos, for 

which the parameters are: Rns = 6.7+5 2 ^ m anc ^ k^eff = 
137tig eV, for M ns = 1-4 M©. However, an addi- 
tional PL component is required to account for the high- 
energy tail of the spectrum, as evidenced by a low f- 
test probability (prob. =0.002). Adding the PL compo- 
nent to nsatmos results in the following improved fit : 
X^/dof (prob.) = 0.97/31 (0.52), with R NS = 6.3+.^ km 
(M NS = 1.4 M Q ), kT cS = 135±geV and V = 2.i±J;| 
(Figure 0). The best- fit projected radius Roo and tem- 
perature are consistent with the qLMXB classification. 
There is no evidence of variability on the timescale of 
the observation as evidenced by a Kolmogorov-Smirnov 



test (probability p = 0.25). 

To improve the statistics, the MOS1 and MOS2 ex- 
tracted spectra are added in XSPEC and a simultane- 
ous spectral fit is performed. A multiplicative constant is 
added to the nsatmos + powerlaw model to verify that 
there is no discrepancy between the pn spectrum and 
the MOS1/2 spectra normalizations. The multiplicative 
constant is fixed to c = 1 for pn and fit for the MOS1/2 
data. It is found that cmosi = l-08+.o.is an d cmos2 = 
0.86+0 .16' both consistent with unity. Consequently, no 
discrepancy is observed and the simultaneous fit is per- 
formed with the values fixed at c = 1 for all three spectra. 
The best-fit nsatmos + powerlaw parameters arc listed 
in Table [3 The NS parameters and the unabsorbed lu- 
minosity (thermal component Lx.th = 4.2xl0 32 erg s _1 , 
and PL component Lx,pl = 2.1 x 10 32 erg s _1 , in the 0.5- 
10 keV range) are within the expected range of values for 
qLMXBs. 

The PL component accounts for 33% of the total unab- 
sorbed flux. Such a strong PL component, if not intrinsic 
to the source, could suggest that, as observ ed for the can- 
didate qLMXB in the core of NGC 6304 (IGuillot et al.l 
2009b), multiple X -ray sources are unresolved in the 
core. Running the detection algorithm wavdetect on the 
MOS1/2 images (because of the smaller pixel size than 
the pn camera) does not reveal the presence of multi- 
ple resolved sources in the core. Nonetheless, visual in- 
spection of the EPIC images seems to suggest that the 
spatial count distribution deviates from that of a point 
source, favoring the hypothesis that nearby sources are 
unresolved by XMM-Newton. 

A more quantitative analysis is performed using the 
XMMSAS task eradial. By comparing the spatial dis- 
tribution of counts to the telescope PSF (FigureEJ), there 
is marginal evidence that the source differs from a point 
source (x 2 =51.5 for 31 degrees of freedom (dof), with 
n.h.p = 1.2%). In other words, there is marginal support 
to the hypothesis that the qLMXB candidate is confused 
with nearby sources (at ~ 15"-30") in the XMM data. 
For this fit and for Figure^ the radial extent of the PSF 
is cut off at 100", the distance of the next closest X-ray 
source detected. However, the fit statistics are highly de- 
pendent on the radial extent of the PSF. For example, 
choosing instead a cut off at 30" leads to x 2 =30.1 for 9 
dof, with n.h.p = 0.04%, which is strong evidence for the 
presence of nearby unresolved sources. 

The hypothesis that the PL component is due to 
nearby sources is further tested by using a spectrum cre- 
ated from counts located within 12" of the source cen- 
ter, therefore decreasing the contamination from possi- 
ble nearby sources at > 12", as evidenced by the PSF 
analysis above. However, this method reduces the S/N 
since only 62% of the ECF is included within 12" at 
1.5 keV. This spectrum has a count rate of 15±1 ctsks -1 . 
Using the nsatmos model alone, the best-fit param- 
eters, R NS = 6.4tJ;gkm (for M NS = 1.4 Mq) and 

kT e s — 135+27 are consistent with the best-fit pa- 
rameters mentioned above, as is the unabsorbed flux 
F x = 1.49xl0~ 13 erg cm -2 s _1 . The statistically ac- 
ceptable fit, xf,/dof (prob.) = 1.73/12 (0.05), does not 
require a PL component, as demonstrated by the high F- 
test probability (prob. = 0.31). Adding a PL component 
leads to consistent NS parameters and flux. Also, the 
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Fig. 3. — XMM surface brightness as a function of radius of the 
candidate qLMXB in the core of NGC 6553. The line indicates 
the best-fit pn-camera point-spread function at the position of the 
source. Counts in excess of the telescope PSF at ~ 15"— 30" are 
evidence of the presence of one (or more) unresolved sources. This 
is marginally supported by the statistic of the PSF fit: x 2 = 51.5 
for 31 d.o.f., with n.h.p.=1.2%. 

fraction of the PL to the total flux is consistent with zero 
at the lcr level. This provides additional support to the 
claim that the PL component represents the emission of 
nearby contaminating sources at angular distances > 12" 
from the qLMXB candidate in the core of NGC 6553. 

Finally, the sec ond XMM-Newton serendipitous source 
catalog f2XMMi: IWatson et all2009f ) reports two sources 
located 20'/7 apart, at the positions R.A.=18 h 09 m 16.5 s 
and decl. = — 25°54'26'.'l for the brightest one, coincident 
with the qLMXB and, R.A.=18 h 09 m 18.0 s and a fainter 
nearby source at decl. = — 25°54'27'/5. The discrepancy 
between 2XMMi and the detections in this work may re- 
side in the algorithm used. While we performed the de- 
tection with wavdetect, the 2XMMi sources are detected 
with a sliding box algorithm. In addition, we used the 
pn data only, while the 2XMMi catalog is constructed 
with a pn+MOSl/2 mosaic. 

Overall, the above analysis supports the hypothesis 
that one or more nearby sources produce the observed 
PL component. However, only the angular resolution of 
the Chandra X-ray Observatory offers the possibility to 
ascertain that the above claim is correct, by spatially 
and spectrally differentiating multiple sources. Archived 
Chandra data of NGC 6553 are analyzed and the results 
are presented in the following subsection. 

3.4.2. Chandra observations of NGC 6553 

The following pertains solely to the confirmation of 
the candidate qLMXB and makes use of Chandra's an- 
gular resolution to understand the core-source emission. 
A systematic detection and spectral analysis of all X-ray 
sources are not reported here. 

Only one source is detected with wavdetect in 
the core of the GC, with a background subtracted 
count number of 46.8±6.9 cts. The position of 
CXOU 180916.48-255426.58 is consistent with the cor- 
rected XMM position of the candidate qLMXB. The dif- 
ference between the Chandra and XMM positions, l'/4, 
is within the accuracy of the Chandra systematic uncer- 
tainty (0'.'6) and the XMM residual error after the cor- 
rection (l'.'l, see Section [2T4]) . Therefore, the Chandra 
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Fig. 4. — XMM contours plotted over the Chandra ACIS-S3 im- 
age. The small circle is the candidate qLMXB detected with Chan- 
dra. The large dashed circle represents the 0.'55 core radius of the 
GC. The red crosses 'x' show the positions of the two 2XMMi 
sources. Finally, the five black '+' signs show the possible sources 
that have low significance in this short ~ 5 ks Chandra observa- 
tions. 

source is consistent with being the same X-ray source as 
XMMU J180916-255425. Figure H shows the Chandra 
image of the GC core, with XMM contours and the two 
2XMMi catalog sources. 

A second possible source, with a 2. lcr significance 
(1.28±0.54ctsks _1 ), lies 17" from the core candidate 
qLMXB, and close (~ 5") to the faint 2XMMi source. 
An additional one, (~ 4") away from the faint 2XMMi 
source and 19'.'5 distant from the qLMXB candidate, has 
a significance of 1.6er. Other faint possible sources, with 
low-significance (1.6c, 1.7c and 2.2c) are located at dis- 
tances 14'.'8, 17'/7 and 22'.'7, respectively (see Figure HJ). 
The extended emission observed in Figure [3] could be 
explained by the presence of other sources in the GC 
core. Unfortunately, this short Chandra observation docs 
not permit a confident detection of all the sources in the 
core. Furthermore, restricting counts to the 3-10 keV or 
the 2-10 keV energy ranges in wavdetect does not lead 
to the detection of any nearby source, suggesting that 
this excess of counts above 2 keV is not due to a sin- 
gle source, but due to multiple fainter sources. A longer 
high angular resolution observation of the central region 
of NGC 6553 would certainly shed light on the PL tail 
observed in the XMM spectrum of the candidate qLMXB 

The PL component measured with XMM has an ab- 
sorbed flux Fx — 4.9xl0~ 14 erg cm -2 s _1 in the 0.5- 
lOkeV range with a photon index T = 2.1. Using 
webPIMMs 7 , the expected Chandra/ ACIS-S3 count rate 
is 4.7 cts ks -1 (0.5-8 keV), which corresponds to 23.5 
counts on this short Chandra observation. However, the 
Chandra background subtracted count rate of an annu- 
lus (r in = 10" and r out = 30") centered around the core 
source is 7 ± 14 cts. The non-detection of an excess of 
counts around the core source does not support the pos- 

Available via HEASARC, http://heasarc.gsfc . nasa . gov/Tools/w3pimms . h 
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TABLE 3 

Spectral analysis of the core candidate qLMXB 



Detector 


Model 




kT afi 


r 


Fx 


PL Flux xl/d-o.f. (prob.) 






(km) 


(eV) 






(96) 



EPIC-pn (12") NSATMOS + PL 6.4±?;J 131±^ 2 -0±i;l l-6li| 23 ±23 1.37/10 (0.18) 

EPIC-pn+MOSl/2 NSATMOS + PL 6.3± ;| 136±|J 2.1±g;| 1.5±£? 33±|| 0.90/54 (0.68) 

EPIC+ACIS-S3 a NSATMOS + PL 6.41^ 1341^ 2.1±g f 1.5±q ? 32±|| 1.02/58 (0.44) 

ACIS-S3 NSATMOS 6. 61^3 134±|| l-Otas cstat=146.98 



Note. — In all the fits, the galactic absorption, modeled with wabs, is represented by a fixed value of 7Vh,22 — 
0.35. r represents the photon index of the power-law (PL) model. .Fx is the total unabsorbed flux, in units of 
10 -13 erg cm" 2 s -1 , in the range 0.5— lOkcV, with ler-crrors. PL Flux (°/ D ) represents the PL contribution to the total 
flux (lrr-errors). 

a For the ACIS-S3 spectrum, the powerlaw normalization is set to zero. The multiplication constant, set free for 
ACIS-S3, is caois = 0.85±° o f e . 



sibility that multiple diffuse sources are responsible for 
the high-energy contamination of the XMM spectrum of 
the candidate qLMXB. This could be explained if the 
sources are variable. 

Alternatively, the detected PL component could be in- 
trinsic to the candidate qLMXB, but display some vari- 
ability, as has been observed for other LMXBs. This 
possibility is investigated by extracting counts in the 
2-10 keV range. With the pn camera, a count rate 
of 0.97±0.31 ctsks -1 is detected, which corresponds to 
0.36 ctsks -1 expected on the ACIS-S3 detector using 
webPIMMs, assuming a photon index L = 2.1. We detect 
1.13±0.47ctsks _1 on the Chandra data in the 2-10 keV 
range. These count rates are consistent and therefore 
there is no observed variability of the high-energy com- 
ponent (> 2keV), at the < 2a level. 

The Chandra spectrum of the candidate qLMXB in the 
core of NGC 6553, i.e., image on-axis, is extracted from 
a 2'.'5 radius circular region to include more than 95% of 
the source energy 8 . For the purpose of the count rate 
calculation only, the background is extracted using a cir- 
cular region of 25" radius around the source , excluding 
5" around the qLMXB itself, which ensures that > 99% 
of the ECF is excluded. 

Due to the small numbers of counts, a binning of 10 
counts per bin is applied for the candidate qLMXB, 
which only guarantees marginally Gaussian uncertainty 
in each bin. The binned spectrum is used simultaneously 
with the three EPIC spectra to evaluate the consistency 
of the spectra. All parameters are tied between the four 
spectra. Only the powerlaw norm is set to zero for the 
Chandra observation. Also, a multiplicative constant is 
added to the model to quantify the variation between 
EPIC and ACIS-S3 spectra. We find that the best-fit 
parameters of the four spectra are consistent with that 
of the EPIC spectra alone (see Tableland Figure [5]) • In 
addition, the multiplication constant, which was left free 
for the ACIS-S3 spectrum, is statistically consistent with 

unity, cacis = 0-85^o 26- This is evidence that, within 
the statistics of the observation, the Chandra spectrum 
is consistent with the three EPIC spectra. 

Finally, we also perform a fit of the Chandra un binned 
spectrum alone, using Cash-statistic (|Cashl [1979T) . The 
background, representing 2.5% of the extracted counts 

8 Chandra Observatory Proposer Guide, Chap. 6, vll.0, Jan 
2009 




0.5 1 2 5 10 

Energy (keV) 



Fig. 5.— XMM-pn (black), MOS1 (red), MOS2 (green) and 
ACIS-S3 (blue) spectra of the candidate qLMXB of NGC 6553, 
XMMU J180916-255425, with the best-fit nsatmos+powerlaw 
model with wabs galactic absorption Nn 22=0-35. The ACIS-S3 
counts are grouped with a minimum of 10 cts per bin, due to 
the small number of counts. The fit statistic x?/dof (prob.) = 
1.02/58 (0.44) and the multiplicative constant applied to the ACIS- 
S3 spectra c = OSbtoil (consistent with unity) are evidence that 
the Chandra spectrum is consistent with the XMM-EPIC spectra, 
within the uncertainty. 

(1.2 counts out of 48 in the source extraction region) can 
be neglected, as required by Cash-statistic. For compar- 
ison purposes with the XMM spectral fitting, the chosen 
model here is nsatmos alone (with fixed wabs absorp- 
tion, distance and mass as before), which assumes that 
the PL component observed in the XMM spectrum is 
not intrinsic to the candidate qLMXB, i.e., that in the 
XMM spectrum, the PL tail is due to contaminating 
nearby sources. The best-fit value are i?NS = 6.6j^' 3 km 
(Af NS = 1.4 M Q ) and kT cS = I34+39CV, which are in 
remarkable agreement with the values obtained from the 
XMM observation. The 0.5-10 keV unabsorbed X-ray 
flux is Fx = 1.0xl0~ 13 erg cm~ 2 s _1 . These results are 
also listed in Table [31 for a convenient comparison with 
the XMM spectral fits. 

4. DISCUSSION AND CONCLUSION 

Using the XMM observation of NGC 6553, we have 
spectrally identified one candidate qLMXB in the GC. 
The best-fit nsatmos parameters found are i?NS — 
6.3tojkm (for A/ NS = 1.4 M ) and kT cS = 136tHeV, 
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and the best-fit PL photon index is T = 2.1±g;|. The 
low-S/N Chandra spectrum finds consistent best-fit val- 
ues within the statistic permitted by the detected counts. 
The modest precision of the NS radius measurement ob- 
tained with these short exposures does not permit to 
place constraints on the dense matter EoS. Deeper Chan- 
dra exposures will provide the S/N necessary to con- 
firm the classification and to produce useful uncertainty 
on the radius measurement for the determination of the 
dense matter EoS. 

The number of candidates identified in NGC 6553 (only 
1) is in agreement with that predicted from the rela- 
tions A qNS ~ 0.04 x r GC + 0.2 dGendre est al.ll2003D and 
N = 0.993 x r GC - 0.046 (jHeinke et al.ll2003bl ). where 
the authors used different normalizations for Toe! the 
encounter rate of the GC is defined by Tqc oc /oJ' 5 r^ 
(jVerbuntJ 12003) . In this last equation, po is the central 
luminosity density (in LqPC -3 ) and r c is the core ra- 
dius (in pc, not angular distance). Quantitatively, us- 
ing the same respective normalizations, the number of 
expected qLMXBs in NGC 6553, given po,ngc 6553 = 
3.84 £opc~ 3 (|Harrislll996t ), are 0.43 and 0.34 using the 
two relations cited above, respectively. Assuming Pois- 
son statistics, the probability of finding one qLMXB 
when 0.43 (0.34) are expected is 28% (24%). Therefore, 
our findings are not in disagreement with these two em- 
pirical predictions. 

The total unabsorbed XMM flux of the candidate 
qLMXB, F x = 1.5±q^x 10~ 13 erg cm" 2 s" 1 , includes a 

PL component with a 33^22% contribution to the flux. If 
intrinsic to the source, this PL component has one of the 
strongest contribution to the total flux ever observed for 
a qLMXB in a GC (others are in NGC 6440 and 47Tuc, 
iHeinke et al.l 12003b. 2005!, respectively). However, we 
have shown evidences for the presence of nearby contam- 
inating sources. This is tentatively supported by a short 
archived Chandra observation, but the low-significance 



detection of the nearby sources does not confirm with 
certainty the existence of contamination. Only a longer 
Chandra observation would confirm the spectral classifi- 
cation of the bright core source and confirm the detection 
of the faint neighboring core sources. 

The candidate qLMXB did not display X-ray variabil- 
ity on the time scale of the observation. Specifically, 
there was no evidence that the XMM integrated light 
curve is significantly different from a linear distribution 
of counts. Regarding the long time scale variability, there 
is no significant change of flux between the XMM and 
Chandra observations. We also searched for archived 
ROSAT observation, but no source was detected at the 
position of XMMU J180916-255425 on the PSPC sur- 
vey, and no HPJ observation at the center of the GC was 
performed. 

To conclude, the candidate in the core of NGC 6553, 
with £ x ,th = 4.2xl0 32 erg s _1 (d = 6kpc), adds to 
the small list of known GC qLMXB. However, only a 
high enough S /N observation with Chandra will allow to 
confirm with certitude the spectral classification of the 
source, and will be able to provide constraints on the 
dense matter EoS. 
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